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We calculate the mass shifts of heavy-light scalar mesons due to hadronic loops 
under the assumption that these vanish for the groundstate heavy- light mesons. The 
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r~| , results show that the masses calculated in quark models can be reduced significantly, 

ji^, We stress that the mass alone is not a signal for a molecular interpretation. Both 

pH ' the resulting mass and the width suggest the observed Dq state could be a dressed 

cq state. We give further predictions for the bottom scalar mesons which can be 
^ [ used to test the dressing mechanism. 
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I. INTRODUCTION 

The constituent quark model has been very successful in describing hadron spectroscopy. 
In recent years, some newly observed hadrons attracted much interest from both the ex- 
perimental and theoretical community since these hadrons do not fit to the quark model 
predictions For instance, the mass of the observed charm-strange scalar meson D^q ^ is 
2317.3±0.6 MeV 13|, while the predictions from most quark models spread from 2400 MeV to 

nn 

2500 MeV [4,15|]. The high mass predicted in constituent quark models is obtained through 
an orbital angular momentum excitation. The result from QCD sum rules in heavy quark 
effective theory gives a mass range 2.42 ±0.13 GeV for the D*q (cs) state which is consistent 
with, but the central value is 100 MeV higher than, the experimental value Q]. Due to the 
fact that the mass of the D*J2317) is iust below the DK threshold at 2.36 GeV, a DK 
molecular interpretation was proposed by Barnes et al. [7| and some others [8|. We want 
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to remark that a DK bound state can be dynamically generated with a mass consistent 
with the observed mass of the i5*g(2317) in the framework of the heavy chiral unitary ap- 



proach 0, [lo|. Other exotic explanations were also proposed, such as tetraquark state [llj . 
and Dtt atom 121]. Besides these exotic explanations, some authors tried to modify the quark 
model predictions. In 13|], one loop corrections to the spin-dependent one-gluon exchange 
potential was considered, and the predicted mass of the D*q is higher than the experimental 
value by only about 20 MeV. Another kind of modification is the mixing of the cs with 
the cqsq tetraquark [l^, or considering the coupling of the cs to hadronic channels, such 
as DK The prediction for the D*q from the QCD sum rules can also be lowered to 
2.331 ±0.016 GeV considering the DK continuum explicitly [l^. Note that it is possible to 



distinguish a hadronic bound state from elementary hadrons, as done for the deuteron 17l | 
and for the light scalars ao(980) and /o(980) [3]. All the results from Refs. [isl, [l6| indicate 
the importance of the strongly coupled hadronic channels on determining the mass of a 
hadron. For the charm-non-strange sector, both the Belle and FOCUS collaborations re- 
ported a scalar meson with a large width |l^, H] . Although the reported masses by different 
collaborations are not consistent with each other, the measurements are considered as the 
same charm scalar meson by the Particle Data Group (PDG), and the PDG average value 
of the mass is 2352 ± 50 MeV. The structure of this state has not been clear yet. In this 
Letter, we revisit the mass shifts of the heavy scalar mesons induced by the strongly coupled 
hadronic loops. For instance, the D^q , the I^Pq cs state in quark model, can couple to the 
D^K^, D^K^ and DsT] loops, see Fig. [TJ The coupling constants of the D^q to the three 
channels can be related by SU(3) symmetry. We shall use three different coupling types to 
study the mass shifts, called Model I, II and III in the following. In Model I, the coupling 
of the scalar heavy meson (S) to the heavy pseudoscalar meson (P) and the Goldstone bo- 
son (0) is assumed to be a constant. In Model II, the coupling is derived in the framework 
of heavy meson chiral perturbation theory (HMxPT) which combines the chiral expansion 
with the heavy quark expansion 2l|, |22] (for a review, see Ref. 23|]). In Model III, a chiral 
effective coupling is constructed disregarding the heavy quark expansion. Of course, such 
corrections due to hadronic loops are always model-dependent and are, in general, difficult 



to quantify, for a recent discussion see 



241 ] (For further discussion of incorporating hadronic 



loops in the quark model, see 25|]). This is why we consider three different models and also 



need to assume that for the ground state meson Qq (with Q = c,b and q = u,d, s), the shift 
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due to the hadronic loops vanishes as it was done e.g. in the calculation of the mass shifts 
of charmonia in Ref. 261 ]. 
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FIG. 1: The relevant hadronic loops coupled to the D*^. 



II. CHOICES OF COUPLING 



Model I 



First, the coupling of S to P,(p is taken as a constant. The loop that integral enters the 
dressed propagator is 



GHs) 



{2n)^ [q^ — m\ + ie)[{p — qY — m\ + ie] 



(1) 



where s = p^. The analytic expression is given by 27|, |28 | 
1 



i? - 1 + In ^ + 



m — 2 + TT I- '■'^ - + m2 + a) 



2s 



167r2 /i2 2s 

— ln(— s + ml — ml + a) + ln(s + ml — ml + a) — ln(— s — ml + m\ + a)] } ,(2) 



where R = — [2/(4 — rf) — 7^; + ln(47r) + 1] will be set to zero in the calculations 
and '-fE is Euler's constant, /i is the scale of dimensional regularization, and a = 
^y[s — (mi + 7712)^] [s — (mi — 7712)^]. Taking into account the D^K^ , D^K^ and D^t] chan- 
nels, the shifted mass of the D*q is given by the solution of the equation 
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(3) 



where M denotes the bare mass of the D*q. Note that we only consider the lowest possible 
intermediate states. In principle, all states with quantum numbers allowed by conservation 
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laws can contribute 



25| . For instance, besides the channels considered here, the D*K*, 



Dgt]' and D*p can contribute either. But their threshold are at least 500 MeV higher than 
that of the DK, thus their contributions are expected to be suppressed. The corresponding 
equation for the non-strange charm meson Dq is 



s-(Mz,2)'-/Re 



6 



G] 



0, 



(4) 



where four channels, D'^tt^, D^tt'^, D^t] and D^K^, are taken into account. The coupling 



constant g has been calculated by using light-cone QCD sum rules in |29|, g = 6.3 ± 1.2 GeV 
in the charm sector, g = 21 ±7 GeV in the bottom sector. A more recent analysis considering 
the Dtn(23n) state as a conventional cs meson gives the coup ling constant for D*qDK as 



GeV 



30|, which is consistent with that given in Ref. 29|]. In this Letter, we study the 



mass shifts of bare cq (and hq) mesons induced by hadronic loops. The values of coupling 



constants given in Ref. 



291 ] will be taken because the masses of the scalar heavy mesons 



used therein are consistent with the quark-model expectation (no fitting to the mass of the 
D*o(2317) was performed since the state had not been discovered yet) and hence correspond 
to the bare masses. 



B. Model II 



In the HMxPT, the Lagrangian for the coupling of S* to P and to leading order is 22|, |23| 



iy/2h 



where the axial field 



contains the Goldstone bosons 



* a=l 



1 J_7rO I J_„ 


7r+ 




















-73 V 



(5) 



(6) 



(7) 



the subscripts a, b represent to the light quark flavor u,d,s, and /^^ = 92.4 MeV is the pion 
decay constant. 



1- 



(P77m - Pal5) 
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represents the multiplet containing the pseudoscalar charm mesons, P = {D^, D^, Df), and 
vector charm mesons, P* = and H = '^qW'^q. 

Sa=^A^ {Dl^.l, - Doa) , (9) 

represents the multiplet containing the scalar charm mesons, Dq = {Dq^ , Dq~^ , D*q) , and 
axial charm mesons, Di = {D^, , D^^). These field operators in Eqs. ([8]) and iQ have 
dimension 3/2 since they contain a factor v^M, where M is the mass of the corresponding 
meson, in their definition. Let p denotes the momentum of a pseudoscalar charm meson. 



e.g. D, and k = p — M^v its residual momentum. The propagator of D in HMxPT is 2l[ | 



where A = M^* — M^- The propagator of the strange charm meson Dg is 



(10) 



(11) 



2{vk + lA.,-6y 

where = M^* — Md and 6 = Mo^ — Md. In the actual calculations, we take 5 = 0.1 GeV 
which is an approximate value of M^^ — Mjj and Mjj* — M^,*. The propagator of the 
scalar charm mesons are similar with proper mass differences, i.e. = M/)^ — M^* for 
Dq, Ass = ^Dsi ~ ^D*o for D*Q, and the SU(3) breaking mass difference can be taken as 
6s = 0.1 GeV, the same as 6. The coupling of the Dq^ to the and 7r~^ to leading order is 



z{7r^q)D'^{q')\C\D*+{p)) = -^^\| M^M o^^v ■ q 



\2 ]\/t2 I ^2 



The last equality holds for on-shell and vr"*" mesons. The relation between h and g is 

h = . ^g. (13) 

Corresponding to g = 6.3±1.2 GeV, calculations via QCD sum rules give h = — 0.44±0.09 for 



charm mesons 



29|. In Ref. 



291] , the mass difference between the scalar and the pseudoscalar 



heavy meson is taken to be 500 MeV. For bottom mesons, we have h = —0.52 ± 0.18 



29|. 



The momentum of a scalar charm meson is p = Msv + k, where Ms is the bare mass of 
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the scalar charm meson and k is the residual momentum. The momentum of the Goldstone 
boson in the loop is denoted by q. Then the residual momentum of the pseudoscalar charm 

o 

meson, whose mass is Mp, in the loop should he k' = p — q — Mpv = k + {Ms — Mp)v — q. 
The loop integral is 

i f d^q {y ■ qY 



G'\v ■ k) 



(14) 



(27^)^ (g2 _ m2 + ie)[{v . {k - q) + M s - Mp + ^ + is] 
where m is the mass of the Goldstone boson in the loop. The loop integral can be worked 



out as 



31 



32| 



G^\v ■ k) 



m 
IGtt' 



■J(0;w) 



(15) 



and 



J(0;cj) = uo{R + \n 



m 



2\/uj'^ — m? cosh ^(— ) — 2'j\%\J uJ^ — m^, u > m 



+ < 



2Vm2 -Cj2cOS^l(-^), 

-2Vio'^ - m2 cosh"^(-^l 



uj < —m 



(16) 



where u = v ■ k + Ms — Mp + 3 A/4. The dressed propagator for D^^ is 



(17) 



2{vk + lAss - Ss) - '-^Mnt^Re [2MnG'iAv " k) + iMn^G^Jv ■ k)] 
For the mass difference A^^ the physical values is taken in the propagator which is correct 

o 

to the order considered. Thus substituting v ■ k hj v ■ p — Mp,*^, the shifted mass can be 
given the value of v ■ p which is the solution of the following equation 311] 



2/^2 o 

2{vp- Mp* 1 - ^Mp* Re 



2MoG'h^{v ■ k) + -Mpfi\^{v . k) 



0. 



The corresponding propagator of the is 



2{vk + \As) - ^Mz^sRe [\MpGl^{v ■ k) + \MpG'h,{v ■ k) + Mpfil^^{v ■ k 



(19) 



From Eq. (ITSl) . the mass shift vanishes in the chiral limit, similar to the mass shift of the 



nucleon due to Nti loop in heavy baryon xPT 
see below. 



3l|, in contrast to that in Model I and III, 
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C. Model III 

Disregarding the heavy quark expansion, we can directly construct an effective chiral 
Lagrangian describing the couphng of a scalar charm meson with a pseudoscalar charm 
meson and a Goldstone boson. The Lagrangian is 

C=h'D,,Ald^Pl + \i.c. (20) 

Note that the field operators of the heavy mesons in Eq. ( l20l) have dimension 1, different 
from those in Eqs. fl8ll9|l . This Lagrangian drives the coupling to be of the type 

t{7.^{q)D\<^)\C\Dl\p)) = ^^q. q' 

. h' 



{MD^Y-Mf,-mi\. (21) 



2^2/. 

The second equality is fulfilled only for on-shell and vr"'" mesons. The relation between 
h' and g is 

h' = '^-^ g. (22) 



For extracting the value of h' , we take M^* — M^i = 500 MeV following Ref. 29|]. Then we 
have h' = 0.78 ±0.15 for charm mesons. Similarly, the coupling constant for bottom mesons 
can be obtained as h' = 1.00 ± 0.33. The dressed propagators of the D^q and Dp~ are 

(23) 



and 

-o ; ' (24) 

s - {Mn*y - ijRe [lG%i{s) + (.) + Ggi^(.)] 
respectively. The physical masses of the charm scalar mesons can be obtained by setting the 
denominators of the propagators to zero. The loop integral in the dressed propagators is 

G^n(s) = ^ f ^ [q-jp-Q)? .25) 

J {271^ {q^ - mj + ie)[{p - qy - ml + le]' ^ ^ 
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The analytic expression can be worked out as 



4 22 4^/2 2\ ^ 



i? + ln^ 



(m^ + m2 — s)^ {mj — 1711 + s) [^^ + 2m^(m^ + m^)] — 2m^cr^ ^ mj 
4 8s m2 

, ^("^1 + "-^i - r, / 2 , 2 , N w ,2 2 , N 
H ^ ym[s — nil + ^2 + o") — ln(— s + — + o") 

+ ln(s + ml~ml + a) - ln(-s - + m2 + cr)] } . (26) 



III. RESULTS 

In general, the effect of a hadronic loop coupling to a bare state is to pull its bare mass 
down to the physical one (if the physical mass is above the threshold of the channel, the 
imaginary part of the loop will contribute to the width of the state). One can expect that 
the mass of a charm meson cannot be pulled down to below the mass of the charm quark 
rric- So we can assume that at some point, called the subtraction point, above rric, the 
contribution from any hadronic loop vanishes. On the other hand, the subtraction point can 
not be too high since a state close to the threshold of a strong decay channel would have 
a strong coupling to the channel, and hence its mass would be affected. A similar idea has 



been taken to study the mass shifts of charmonia [26|]. After choosing a specific subtraction 
point, the renormalization scale /i considered as a parameter can be determined from this 
assumption. In the loop function in Model I, Eq. ([2]), the coefficient of the chiral logarithm 
In (m^/yU^) is a constant l/(167r^); in the loop functions in Model II and III, Eq. f|T5|l and 
Eq. (1261) . the coefficients of the same logarithm are momentum-dependent. In Model III, 
the coefficient 



III 



167r2 



4 22 4*^/2 2\ -^2 

+ m^m2 + ^2 — ^(^1 + ^2)^ + ^■^ 



is always positive, and it changes slowly with respect to ^/s below 3 GeV for the DK loop, 
see Fig. [2](a). So in Model I and III, the dependence of fi on the choice of subtraction point 
is small. For instance, the values of fi determined when the subtraction point is chosen at 
^/s = Mo = 1.87 GeV and ^/s = rric = 1-35 GeV are listed in Table [B One can see when 
the subtraction point is changed from to rric, the resulting values of fi in the DK loop 
changes slightly, and all the values are not far away from Md. In the following, we shall 



9 



take Mf) as the subtraction point for Model I and III. However, in Model II, the coefficient 
Cii = m^Ci;/(167r^) is proportional to v-p, and changes its sign at just below Md. In Fig.[2](b), 
we show Cii as a function of v ■ p where Mp = Md and m = rriK are used. So the resulting 
value of /i depends strongly on the choice of subtraction point. The values of fi in Model II 
are also given in Tabled It seems that these values are unphysical. Because we only take the 
leading order in heavy quark expansion to describe the coupling in Model II, the unphysical 
value of /i might indicate that higher order contributions, which will change the behavior of 
Cii, are important. To avoid this problem, a commonly used value /i = 1 GeV in HM^PT 
is taken in Model II. Certainly, the /i dependence should be absorbed by counterterms at 
the next order. The same method can be used directly for the scalar bottom mesons 
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FIG. 2: The coefficients of the logarithm In (?7T,2//i^) in the loop functions in Model III (a) and 
Model II (b). nil = and m2 = mx are taken. 



TABLE I: The values of the renormalization scale /x in the DK loop determined when the subtrac- 



tion point being chosen at AId and 


rric- All units are 


in GeV. 




Subtraction point 


Model I 


Model II 


Model III 


Md = 1.87 GeV 


1.58 


1095.3 


1.85 


rric = 1.35 GeV 


1.84 


0.20 


1.82 



with no more free parameters. In Model I and III, we use the same assumption that the 
contribution from any hadronic loop vanish at the mass of the lowest heavy flavor meson. 
That is, in the bottom case the subtraction point is chosen to be Mb = 5.28 GeV. The bare 
masses of the scalar heavy mesons are taken from the popular Godfrey-Isgur quark model 
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which can describe the meson spectroscopy, especially the low-lying states very well 
Md' = 2.48 GeV, Md- = 2.40 GeV, Mb* = 5.83 GeV and Mb* = 5.76 GeV. 

The resulting masses of all the four mesons are listed in Table HTl The results from Model I 

TABLE II: The resulted masses of the lowest scalar heavy mesons from dressing in different models. 







Bare Mass 


Model I 


Model II 


Model III 


AMiii (GeV) 




(GeV) 


2.48 


2.33-2.39 


2.26-2.35 


2.36-2.40 


0.08-0.12 




(GeV) 


2.40 


2.30-2.35 


2.39-2.40 


2.30-2.34 


0.06-0.10 


Mb* 


(GeV) 


5.83 


5.58-5.72 


5.73-5.73 


5.62-5.70 


0.13-0.21 


Mb* 


(GeV) 


5.76 


5.49-5.67 


5.85-5.86 


5.55-5.64 


0.12-0.21 



are consistent with those in Model III, which show that the bare mass of the D*^ can be 
pulled down significantly to the region close to the mass of the D*q(2317) state. If we choose 
a bare mass from another quark model, the obtained mass can even be consistent with the 
experimental value. In Model I, a constant coupling is taken which would violate Goldstone's 
theorem because the vr, K and rj are Goldstone bosons. Model II has a large // dependence, 
which makes it not preferable for a phenomenological analysis. The Lagrangian for Model III 
is constructed from chiral symmetry, and the /i dependence is really small (for instance, the 
resulting mass of the -D*q would change to 2.37-2.41 GeV if we choose rric = 1.35 GeV as 
the subtraction point, and the change is no more than 10 MeV). Therefore, we choose the 
results in Model III to give further predictions. The absolute mass shifts in Model III are 
listed in the last column in Table [TTl The mass of the DX is consistent with the experimental 

n 

data 2352 ± 50 MeV [3|]. One strong decay channel Dtt is open for the state, and the decay 
width of this channel should give the dominant contribution to the width of the Dq. Then 
the width of the Dq can be obtained from 

r^. = 5' -ImG"l(M^.)- (27) 

^ If the parameters of a quark model, e.g. the constituent quark masses and strong coupling constant, were 
determined by fitting to the whole hadron mass spectrum, loop effects would be incorporated to some 
extent although in an unclear way. But in reality, the quark model parameters were fitted to mostly the 
low lying states. Especially in the Godfrey-Isgur quark model the physical spectrum used in fit does not 
contain the heavy scalar mesons studied here. 
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Using the mass M^)* from Table [TTl we obtain 

r^* = 99 - 167 MeV . (28) 
The result is roughly consistent with the experimental width for the Dq, which was reported 



as 276 ± 21 ± 63 MeV by the Belle Collaboration [19| and 240 ± 55 ± 59 MeV by the FOCUS 



Collaboration (20|]. Both the mass and the width suggest the observed Dq state can be the 
dressed cq state. Note that in Ref. [l^, two molecular states were predicted, and both of 
them are not consistent with the data. Certainly, the uncertainty of the data is large so 
far, and more precise measurements are highly desirable. Another noticeable result is that 
the mass shifts in the bottom sector are about twice of those in the charm sector, as can 
be seen from the last column in Table [TTl Based on heavy quark symmetry, and assuming 
the D*q(2317) is a DK bound state, the mass of the BK bound state was predicted as 
5733 MeV 33 1 which was confirmed by a dynamical calculation |l^, larger than the mass 



region obtained here by dressing the bq state. That means if a B*q state with a mass which 
is much smaller than 5733 MeV were found, it would probably be a dressed bq state rather 
than a BK bound state, or the i5*Q(2317) would not be a DK bound state. Similar to that 



of the Z)q, the width of the Bq can be estimated as 



62 - 100 MeV . (29) 



IV. SUMMARY 



The D*o(2317) is considered as a DK molecular state by many authors because in the 
quark model its mass simply comes out too high. However, hadronic loops can pull down its 
mass. To make the calculation of such an effect quantitative, we assume that the hadronic- 
loop induced mass shifts of a hadron vanish at some point. The point is chosen as ^/s = m£, 
for the charm sector and irts for the bottom sector. Then we calculate the mass shifts of 
the heavy scalar mesons by using three different types of coupling. The input bare masses 
are taken from the Godfrey-Isgur quark model ^. The mass of the D*q state is lowered 
significantly, and it can even be pulled down to 2317 MeV if we adjust the bare mass in 
the predicted region from different quark models. That means the simple argument against 
the quark model for too high masses is not valid. What the quark model predicts is just 
the bare masses of hadrons. In order to compare with the experimental spectroscopy, the 
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bare hadron masses need to be dressed 



24| . Note, however, that such a dressing is always 



model-dependent and must be considered in different approaches, as done here. Because 
chiral symmetry is fulfilled and the /i dependence is small, we choose Model III to give 
further predictions. The results of both mass and width show that the observed Dq could 
be a dressed cq state. We also give predictions for the dressed bottom scalar mesons. The 
mass of the B*q is smaller than that of the BK bound state which was obtained assuming 
the D*q(2317) be a DK bound state 3^. Precise experimental data from B factories are 
highly desirable to test the dressing mechanism. We would like to stress that a consistent 
treatment of the mass and width is required. The difference from the molecular state will 
presumably be revealed in the decay pattern with into various channels, this is a further 
step to be investigated. 
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